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ABSTRACT 



We present new sensitive wide-band measurements of the fine structure fine 
'^Pi ^Pq (J=1-0, 492 GHz) of neutral atomic carbon (CI) in the two typ- 
ical Ultra Luminous Infrared Galaxies NGC 6240 and Arp 220. We then use 
them along with several other CI measurements in similar objects found in the 
literature to estimate their global molecular gas content under the assumption 
of a full CI-H2 concomitance. We find excellent agreement between the H2 gas 
mass estimated with this method and the standard methods using ^^CO. This 
may provide a new way to measure H2 gas mass in galaxies, and one which may 
be very valuable in ULIRGs since in such systems the bright ^^CO emission is 
known to systematically overestimate the gas mass while their ^^CO emission 
(an often-used alternative) is usually very weak. At redshifts z > 1 the CI J=l-0 
line shifts to much more favorable atmospheric windows and can become a viable 
alternative tracer of the H2 gas fueling starburst events in the distant Universe. 

Subject headings: galaxies: starbursts - galaxies: individual (NGC 6240, Arp 220) 
- ISM: atoms - ISM: molecules 
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1. Introduction 

The CI J=l-0 line emission emerging from molecular clouds has been detected early 
(Phillips & Huggins 1981), but its interpetation as emanating from only a narrow CII/CI/CO 
transition zone at the surface of these clouds (e.g. Tielens & HoUenbach 1985a,b) seems to 
preclude its use as a bulk H2 gas tracer, except perhaps to derive lower limits. Nevertheless, 
despite the considerable difficulties involved in imaging wide areas in CI emission (caused 
by an unfavorable atmospheric window at 492 GHz, and the small beams of single-pixel 
detectors mounted on large aperture sub-mm telescopes), early work has shown CI to be 
widely distributed and well-correlated with other H2 gas mass tracers such as ^^CO (Keene 
1985; Keene 1997 and references therein). Further observational efforts have confirmed this 
CO-CI concomitance in molecular clouds together with a remarkably constant N(C)/N(CI) 
ratio over the bulk of the CO-bright gas (Ojha et al. 2001; Ikeda et al. 2002). A new 
comprehensive study has shown how this could be estabhshed and advocated the use of CI 
hne emission as a molecular gas mass tracer in galaxies (Papadopoulos, Thi, & Viti 2004). 

Ultra Luminous Infrared Galaxies (ULIRGs) are powered mainly by starburst events 
(Genzel et al. 1998) that yield Lfir > lO^^L© making them the most luminous objects in 
the local Universe (Sanders et al. 1988; Sanders & Mirabel 1996). A number of reasons make 
them a good sample to test CI emission as a bulk H2 gas mass tracer in galaxies, namely: 

a) They are gas-rich systems but with their H2 distributions usually so compact (e.g. Downes 
& Solomon 1998; Bryant & Scoville 1999) that single pointings with single-dish telescopes 
are adequate to record all their CI/CO flux density (of prime importance for the difficult CI 
observations) . 

b) Their H2 gas mass has been estimated using several methods other than the standard 
one based solely on the ^^CO J=l-0 line luminosity (e.g. Solomon et al. 1997; Downes & 
Solomon 1998; Glenn & Hunter 2001; Lihong et al. 2003). Thus they constitute an ideal 
test-bed for any new H2 mass-tracing method. 

c) In these merger/starburst systems intense FUV radiation and powerful tidal fields drive a 
two-phase differentation of the molecular gas (Aalto et al. 1995). The diffuse, ^^CO-bright, 
and non-virialized intercloud phase is probably responsible for the systematic overestimate 
of H2 gas mass by factors of ~ 4 — 5 in starburst environments (Solomon et al. 1997; Downes 
& Solomon 1998). It is thus interesting to see how a new method of measuring H2 gas mass 
will fare under such conditions. 

d) Early applications of this method in Arp 220 (Gerin & PhiUips 1998), and the Cloverleaf 
QSO at z— 2.5 (Weiss et al. 2003) have yielded an H2 mass that is in good agreement with 
the standard methods that use ^^CO and dust continuum emission. 
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e) ULIRGs, as possible precursors of optically-bright QSOs (Sanders et al. 1988; Tacconi 
et al. 2002), or as the nearby templates for starburst, dust-enshrouded galaxies observed 
at high redshifts (Small, Ivison & Blain 1997; Hughes et al. 1998), are interesting systems 
in their own right. Molecular gas being the "fuel" of their prodigious star-forming activity 
makes any new method of estimating its mass particularly valuable. 

In this work we present new sensitive observations of the CI J=l-0 line in two prominent 
ULIRGs namely Arp 220 and NGC 6240, and collect available CI measurements in similar 
objects from the literature, which are then used to estimate their H2 gas mass content. 
Throughout this paper we adopt an Ho = 75 km s~^ Mpc~^ and qo = 1/2. 

2. Observations and results 

We used the 15-mctcr James Clerk Maxwel Telescope (JCMT) on Mauna Kea, Hawaii, 
to observe the ^Pi ^^Pq (J=1-0) fine structure line at u^^st = 492.160 GHz in the ULIRGs 
Arp 220 (January 3nd 2003) and NGC 6240 (February 27th 1999). Cold and extremely dry 
weather with T225GHZ ^ 0.035 ensured excellent conditions for high-frequency observations 
with overall system temperatures of the SSB-tuned W-receiver of Tgyg ~ (700 - 900) K. 
Rapid beam-switching at a chop-throw of 30" (Az) with a frequency of 1 Hz produced stable 
basehnes, and the use of the Digital Autocorrelation Spectrometer (DAS) at its widest mode 
of 1.8 GHz (~ 1100 km s~^) enabled the full velocity coverage of the wide extragalactic lines 
often observed in such galaxies. This is particularly important since e.g. Arp 220 has a 
FWHM ~ 500 km s~^ and in the past multiple frequency settings of the spectrometer were 
necessary in order to cover the entire CI J=l-0 line (Gerin & Phillips 1998). The aperture 
efficiency was measured repeatedly using Mars, an important task in sub-mm telescopes 
operating at such high frequencies since residual thermal distortions and mechanical defor- 
mations of the dish can then affect their performance significantly. We deduced (rja) — 0.35, 
lower than the one quoted in JCMT Observer's Guide (77a = 0.45), but within the expected 
variations (Priberg, private communication). The beam size at these frequencies (used to 
estimate the aforementioned efficiencies) is ^hpbw = 10". Frequent pointing checks showed 
the pointing error to be ~ 3.5". 

Frequent observations of strong spectral-line sources yielded a calibration uncertainty 
of ~ 20%, which along with the ~ 25% uncertainty of adopted value of aperture efficiency 
dominate the error in the reported velocity-integrated line fluxes. These are estimated from 

S„=/ S„dv=Ji^K.(x)/ TXdv = 1^^MJ2k„(x) / Tldv, (1) 
Jay ?7a7rD^ Jay Va. Jay 
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where K,{x) = xV(l - e-^'),x = ^s/(1.2^hpbw) {0 s=source diameter) accounts for the 
geometric coupling of the gaussian part of the beam with a finite-sized, disk-hke source. For 
NGC 6240 it is K^{x.) = 1.15, while for all other sources K^{x) - 1. 

The CI spectra are shown in Figures 1, 2 along with overlayed ^^CO and ^^CO J=2-l 
spectral lines. The correspondance is excellent, while the weakness of the observed ^^CO 
emission (an often-used alternative to the much more optically thick ^^CO) elevates CI 
emission (if indeed fully concomitant with H2) to a good alternative H2 gas tracer in ULIRGs. 

3. Global molecular gas mass estimates using the CI line emission 

We use the measured Sci to deduce the global molecular gas mass from 
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where h = 0.75, the Einstein coefficient Aio = 7.93 x lO-^s"^, and Qio = Qio(n, T^) depends 
on the gas excitation conditions (see Papadopoulos et al. 2004). The [CI]/[H2] abundance 
chosen here is Xci = 3 x 10-^, identical to that used by Weiss et al. (2003) in their H2 mass 
estimate in the Cloverleaf. This value is the average between the minimum found for the bulk 
of the Orion A and B clouds (~ 10^^; Ikeda et al. 2002) and in the starburst environment 
of the nucleus of M82 (~ 5 x 10"^; White et al. 1994) (assuming [C0]/[H2] = 10"^). The 
results are shown in Table 1 where Mci(H2) can be directly compared with the best available 
Mco(H2) estimates. In the cases where a range of physical conditions is reported for the 
gas, the average Qio excitation factor is adopted. 

The agreement is remarkable, especially given that the CO-derived H2 masses reported 
here are based on much better methods than the mere application of the Galactic X factor 
(which in ULIRGs overestimates H2 mass by a factor of ~ 5). It must be noted though that 
even such methods leave uncertainties of a factor of ~ 2 in their M(H2) estimates (e.g. see 
Tacconi et al. 1999). In the one case of IRAS 10565+2448 it is Mci(H2)/Mco(H2) ~ 4, 
possibly too large a difference to be attributed solely to the aforementioned uncertainties of 
the CO-based methods. In starburst environments higher Xci values are possibly more ap- 
propriate (e.g. Schilke et al. 1993 and references therein). Indeed, assuming CI as abundant 
as in the center of M 82 would keep the CI and CO-deduced H2 mass estimates in agreement 
within factors of ~ 2, while bringing Mci(H2) in IRAS 10565-1-2448 in better accord to the 
CO-deduced value. Nevertheless the latter case highlights an irreducible uncertainty when 
using optically thin emission from any molecular or atomic species to trace H2 gas mass. 
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It is important to note that unlike quiescent spirals ULIRGs can have > 50% of their H2 
mass at densities n(H2) > 10^ cm~^ (Gao & Solomon 2004), for which one-dimensional static 
PDR models predict much lower Xci values, effectively rendering CI emission incapable of 
tracing its mass. This does not seem to be happening, lending further support to a full 
CI-II2 concomitance over most of the parameter space characterizing the physical conditions 
of H2 in galaxies. Furthermore CI emission may be a good alternative molecular gas tracer 
particularly in ULIRGs because the often-used optically thin ^^CO emission is very weak, 
while ^^CO J=l-0 and a Galactic conversion factor overestimates M(H2) because of the 
presence of non-viriahzed ^^CO-bright H2 gas. Apart from the uncertainty of the assumed 
Xci (which plagues also the ^^CO measurements since C and ^^CO emerge from tightly 
coupled chemical routes), the CI J=l-0 line is a more straightforward II2 mass tracer than 
CO due to a simpler partition function and a significantly reduced sensitivity to the ambient 
gas excitation conditions. 

It can be argued that a proper comparison between the various tracers of molecular gas 
must be between transitions within the same frequency domain because of the significant 
variations in the atmospheric absorption (and thus in the resulting effective system sensitiv- 
ity) across the mm/sub-mm wavelengths. For CI and CO this is presented by Papadopoulos 
et al. (2004) where it is shown that while the ^^CO J=4-3 (461 GHz) transition can be indeed 
much stronger than CI J=l-0, the latter offers decisive advantages for cooler and/or diffuse 
gas (Tk < 20 K, n < 10^ cm~^) where the emission of the former is significantly surpressed. 
In ULIRGs the very weak ^^CO emission makes the much brighter CI J=l-0 a better alter- 
native H2 mass tracer when using the same telescope to observe unresolved sources. This 
of course assumes CI observations conducted during suitable dry weather conditions, which 
occur much less frequently than the those needed for sensitive ^'^CO observations. However 
this situation is expected to change drastically when the new generation of mm/sub-mm 
single dish telescopes (APEX) and interferometer arrays (ALMA) become operational in the 
excellent sites in the Atacama Desert plateau in Northern Chile. There dry atmospheric 
conditions, suitable for CI J=l-0 observations, are expected for ~ 50% of the time and with 
much less pronounced diurnal variations (Radford & Nyman 2001 and references therein). 

More CI J=l-0 observations of ULIRGs are urgently needed to further test the notion 
of CI as a bulk H2 gas mass tracer. Application of this technique to quiescent nearby galaxies 
is also of prime importance but then extensive CI imaging may be needed to match the CO 
and CI imaged areas. Finally CI observations of galaxies at high redshifts with known CO 
detections can further test the effectiveness of this method while circumventing the problem 
of the low atmospheric transmission at 492 GHz (J=l-0) and 809 GHz (J=2-l). 
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4. Conclusions 

We report on new sensitive CI J=l-0 measurements in the two prominent ULIRGs 
NGC 6240 and Arp 220, and collect CI measurements available for similar objects from 
the literature. We then use them to derive the molecular gas mass under the assumption 
of fully concomitant CI, CO, and H2 and find very good agreement between the CI and 
the CO-based estimates within the uncertainties of the two methods. The faintness of the 
^^CO emission (often used instead of the more optically thick ^^CO) in these two galaxies 
and in similar systems makes CI emission a good alternative tracer of their molecular gas 
mass. However given the much stricter constraints on the atmospheric conditions needed for 
sensitive CI observations (very dry weather) than those for ^^CO, the practical advantages of 
CI as a good optically thin II2 gas mass tracer will be fully realised when the new generation 
of mm/sub-mm single-dish (APEX) and interferometer arrays (ALMA) become operational 
in Cerro Chajnantor in the Atacama Desert plateau. 

We advocate more CI ,1=1-0 observations of similar objects at low as well as high 
redshifts (where CI lines are redshifted to more transparent atmospheric windows) in order 
to further test the potential of CI to trace bulk H2 gas mass in galaxies. 

We thank the referee for comments on the original manuscript that greatly improved the 
presentation of the main thesis of this work. TRG acknowledges support from the Danish 
Research Council and from the European Union RTN network, POE 
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Fig. 1.— Arp 220 (ai95o : 15'' 32"' 46.73^ ^1950 : +23° 40' 07.36"). CI J=l-0 (thick line), 
i^co J=2-l (strongest spectral line), 10 x [^^CO J = 2-1] (thin line). The velocity-averaged 
^^CO/^^CO intensity ratio is R21 ~ 35. The estimated rms noise of the CI spectrum is 
6T\ — 10 mK at a velocity resolution of AVchan = 7.6 km s~^. 
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Fig. 2.— NGC 6240 (aigso : 10^ 21'° 31. P, ^1950 : +47° 24' 23.0"). CI J=l-0 (thick line), 
i^co J=2-l (strongest spectral line) 10 x [^^CO J = 2 — 1] (thin line). The velocity-averaged 
^^CO/^^CO intensity ratio is R21 ~ 50. The estimated rms noise of the CI spectrum is 
5Ta = 8 mK at a velocity resolution of AVchan = 7.6 km s""^. 
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Table 1. Molecular gas mass estimates based on CI and CO 



Name 


Redshift 


Sci 


(Tk,n)b 


Qio 


Mci(H2) 


Mco(H2)^ 






(Jy km s~^) 


(K, cm-3) 




(xlO'^ M0) 


(xlO^ M0) 


NGC 6240^^ 


0.024 


600±180 


(59, 10^) 


0.52 


(6 ±2) 


4 


Arp 220 


0.018 


1160±350 


(64, 10^-10^) 


0.47 


(7 ±2) 


5 


Mrk 231 


0.042 


200±55^ 


(74, 3.5x103) 


0.45 


(7 ±2) 


4 


NGC 6090^ 


0.029 


270±55^ 


(40, 10^-10^) 


0.48 


(4±1) 


6 


10565+2448 


0.043 


500±150^ 


(60, 1.2x10^) 


0.51 


(16 ±5) 


4 



^ from Gerin & Phillips 2000 

^ Solomon et al. 1997; Downes & Solomon 1998 (except NGC 6240, NGC 6090) 
^ Downes & Solomon 1998 (except NGC 6240, NGC 6090) 
Mco(H2) from Tacconi et al. 1999; (T^, n) from Large Velocity Gradient modeling of the 



line ratios reported by Tacconi et al. and references therein, 
c CO J=l-0 flux from Bryant & Scoville 1999; Zhu et al 1999, Mco(H2) deduced by 
using the X factor advocated by Downes & Solomon 1998. 



